Effects of ^He Impurity on Solid ^He Studied by Compound Torsional Oscillator 
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Frequency shifts and dissipations of a compound torsional oscillator induced by solid *He samples 
containing '^He impurity concentrations (x^ = 0.3, 3, 6, 12 and 25 in units of 10~®) have been 
measured at two resonant mode frequencies (/i — 493 and /2 = 1164 Hz) at temperatures (T) 
between 0.02 and 1.1 K. The fractional frequency shifts of the /i mode were much smaller than 
those of the /2 mode. The observed frequency shifts continued to decrease as T was increased 
above 0.3 K, and the conventional non-classical rotation inertia fraction was not well defined in all 
samples with X3 > 3 ppm. Temperatures where peaks in dissipation of the mode occurred were 
higher than those of the /i mode in all samples. The peak dissipation magnitudes of the /i mode 
was greater than those of the /2 mode in all samples. The activation energy and the characteristic 
time (to) were extracted for each sample from an Arrhenius plot between mode frequencies and 
inverse peak temperatures. The average activation energy among all samples was 430 mK, and tq 
ranged from 2xl0~^ s to 5xl0~"^ s in samples with X3 = 0.3 to 25 ppm. The characteristic time 
increased in proportion to x^^^ . Observed temperature dependence of dissipation were consistent 
with those expected from a simple Debye relajcation model if the dissipation peak magnitude was 
separately adjusted for each mode. Observed frequency shifts were greater than those expected from 
the model. The discrepancies between the observed and the model frequency shifts increased at the 
higher frequency mode. 

PACS numbers: 67.80.-s 
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I. INTRODUCTION 

The discoveryi that the resonant frequencies of tor- 
sional oscillators (TOs) containing solid ''He samples in- 
creased at temperatures below about 300 mK has stirred 
much excitement (see reviews^^— ). Interpreted as a par- 
tial decoupling or so-called^ "non-classical rotational in- 
ertia (NCRI)" of the solid ^He samples, the discovery 
gave evidence for the long-sought^ supersolid state of 
solid ''He possessing simultaneously crystallinity and su- 
perfluid properties. In the temperature range where the 
observed TO frequency shifts vary most rapidly, peaks in 
the energy dissipations of the TO are also observed. The 
frequency shifts and the accompanying dissipations of 
TO have been confirmed in many laboratories. ^"^'^ How- 
ever, the interpretation in terms of supersolid state in 
the loaded solid "^He samples remains controversial and a 
comprehensive understanding of the TO frequency shifts 
and other related phenomena in quantum solid ^He at 
low temperatures has not been established. 

There are puzzling observations in experiments de- 
signed to test the supersolid interpretation. Superflows 
that would be expected to occur through supersolid ^He 
samples under applied pressure gradients have not been 
observed in the earlier flow experimentsA^^il carried out 
in the same temperature range as the NCRI effect is ob- 
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served. Recently, however, unusual mass flow; 
induced through solid *He samples below about 0.6 K by 
applying chemical potential gradients across the samples. 
Propagation of fourth sound that would be expecteS^^ 
for a superfluid has not been detectedJ^i^ The shear 
modulus of a thin solid ''He slab was found^^ to "stiffen" 
in a very similar manner as the frequency of TO contain- 
ing an annular solid *He sample was found^ to increase 



with decreasing temperature. This similarity indicated a 
common origin of these two effects. While almost none of 
these observations are clearly understood, it has become 
apparent that the details of the observations are affected 
strongly by the solid ^He sample quality depending on 
growth condition, sample geometry and size, and '^He 
impurity. The objective of the present work is to gain 
understanding of the puzzling role played by '^He impu- 
rity by use of our compound TO techniques. 

The surprisingly high sensitivity of observed frequency 
shifts to minute ^He impurity concentrations (x^) at 
parts per million (ppm) levels in the ^He samples has 
been reported by Kim, et aljii^i The temperature depen- 
dence of the shear modulus was also foundl^ to depend 
on ^3. On the other hand, a thermodynamic anomaly 
in heat capacity occurs at temperatures which are rel- 
atively insensitive to the value of X3. The dynamics of 
■^He atoms within solid "^He may be probed by NMR ex- 
periments. Simultaneous observations of NMR and TO 
effects, of the same crystal, have been carried out by 
Toda, et ali^i^ in solid ^He samples with X3 down to 10 
ppm. They found three different spin-lattice relaxation 
times suggesting the existence of three different states of 
■^He atoms in the solid matrix of "^He. Kim, et al.— 1^ 
reported their measurements of NMR relaxation times, 
Ti and T2, of ^He contained in solid ^He samples with 
X3 down to 16 ppm. The relation between all of these 
NMR observations and the role of "^He impurity in the 
supersolid phenomenon is yet to be clarified. 

One goal of our work is to measure how the dissipations 
accompanying frequency shifts of solid-^ He-loaded TO 
depend on the amount of added "^Hc impurity. The earlier 
reportaii^i focused mainly on the variation of frequency 
shifts as X3 was changed. Another goal is to measure 



how the dependence of frequency shifts and dissipations 
on the TO frequency varies as is changed. Frequency 
dependence studies are difficult in single resonance mode 
TOs. We report on the first systematic measurements 
of the dependence of dissipations and frequency shifts 
on X3 (0.3 ^ 25 ppm) using our compound TO having 
two resonance mode frequencies (~500 and ^1200 Hz). 
Unlike in the earlier ^He impurity dependence study^i^ 
all samples are grown in the same torsional oscillator fol- 
lowing nearly identical solid growth procedures such that 
variance due to sample cell geometry and sample qual- 
ity would be minimized. Measurements are made simul- 
taneously at the two frequencies under identical sample 
conditions and they allow us to explore dynamical ef- 
fects which cannot be probed in a single mode TO. For a 
given sample, the maximum in dissipation of the higher 
frequency mode appears at a higher (peak) temperature 
than the lower frequency mode. As X3 is increased, the 
peak temperatures for both modes increase. Arrhenius 
plots of the mode frequencies vs. inverse peak tempera- 
tures are analyzed to extract the activation energy and 
the characteristic time involved in the dissipation pro- 
cess. A simple form of Debye dissipation combined with 
the extracted characteristic time is inadequate to rep- 
resent the observed frequency dependent dissipations of 
the two modes separately. The frequency shifts expected 
from the dissipation by the Kramers-Kronig relation can- 
not account for those observed in the two modes. The 
limiting characteristic time tq extracted from the two 
modes varies in proportion to x^^^. 



II. EXPERIMENT 

Our compound torsional oscillator shown in Fig. [T] 
was modified from the earlier one^°'^® and has two tor- 
sion rods (a, b) each with diameter 1.9 mm and length 15 
mm, and two interconnected masses. The flange above 
the upper rod (a) is rigidly attached to a large copper 
vibration isolation block in good thermal contact with 
the mixing chamber of a dilution refrigerator. The upper 
mass is a "dummy" comprised of a central disc (c) and 
two electrode fins (d) made of aluminum plates. The TO 
is capacitively driven by applying dc-biased sinusoidal 
voltage between a stationary electrode (not shown) and 
one of the movable electrode fins. The motion of the TO 
is detected by measuring the voltage induced between the 
other electrode fin biased against another stationary elec- 
trode (not shown) . The lower fin (e) acts as an auxiliary 
electrode for measuring the motion of the sample cham- 
ber itself. The lower mass (f) is mostly made of Stycast 
1266 epoxy^^. cast around the base (g) below the lower 
torsion rod. The sample chamber (h) for solid ''He is an 
annular space (8.0 mm inner diameter, 10.0 mm outer 
diameter and 8.0 mm height). Helium is introduced into 
the chamber via the fill hole (i) drilled (diameter = 0.8 
mm) through the center of torsion rods and the base, 
and a diametrical channel just below the lower surface of 
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FIG. 1. Schematic of compound torsional oscillator, (a) upper 
torsion rod, (b) lower torsion rod, (c) upper disc, (d) electrode 
fins, (e)lower electrode fin, (f) cell body, (g) base, (h) annular 
sample chamber, (i) fill hole, (j) copper foil liner. The shaded 
region is machined from a single block of BeCu. The width of 
the annular sample space is 1.0 mm but shown exaggerated 
for clarity. 



the base. To improve the thermal contact between the 
sample helium and the mixing chamber, the lower sur- 
face of the BeCu base is pressed against a 100 /im thick 
copper foil (j) which is extended onto the inner wall of 
the sample space. 

The ■^He impurity concentration is increased by adding 
a calibrated amount of "^He gas into the cell prior to load- 
ing it with the commercial ultra high purity *He (nomi- 
nal X3 = 0.3 ppm) at 4.2 K to approximately 80 bar. A 
solid sample is subsequently grown by the blocked capil- 
lary method. The solid plug formed in the rapidly cooled 
portion of the filling capillary maintains a constant mass 
in the sample space below the plug. A Straty- Adams ca- 
pacitive pressure sensor attached to the isolation block 
is used to monitor the pressure in the fill tube during 
sample solid formation. The measured freezing tempera- 
ture where a sudden increase in the oscillation amplitude 
occurs gives the pressure in the solid formed. After the 
plug is formed, the total time elapsed to freeze ^He in 
the sample chamber is about 40 minutes for all the sam- 
ples reported here. Since the sample solid pressure is 40 
bar or greater, where the freezing temperature is higher 
than the highest superfluid *He transition temperature, 
redistribution of '^He by the heat flush effect is expected 
to be unimportant. Upon completion of measurements 
for a given sample, the dilution refrigerator system is 
warmed up to near 6 K for pumping out helium from 
the cell. The pressure and the TO resonant frequencies 
are monitored during this "bake out" procedure to en- 
sure that the remaining residual amount of helium gas 
is sufficiently small. Subsequent to this procedure, an- 
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other calibrated amount of ^He is added such that the 
■^He impurity concentration is larger than the previous 
sample and the above procedure is repeated to grow the 
next solid sample. In all, measurements were made on 
solid ^He samples with 2:3 = 0.3, 3, 6, 12 and 25 ppm. 
The samples are identified by their '^He impurity concen- 
trations. The uncertainty in X3 is estimated to be ±20 
%. 

Sample temperature (T) is inferred from a ruthenium 
oxide resistance thermometer attached to the isolation 
block and is calibrated against a "^Hc melting pressure 
thermometer. Reproducibility of the ruthenium oxide 
thermometer upon cycling up to room temperature is ver- 
ified with a fixed point superconducting standard down 
to 15 mK. 

The two modes of the compound TO may be excited 
and detected simultaneously. The resonant frequency 
and the amplitude of each mode are tracked continuously 
and independently by two automatic phase-lock feedback 
data acquisition systems. In all of the data presented in 
this report, measurements of both modes are taken si- 
multaneously to ensure identical sample conditions. The 
drive levels are set such that sample velocity amplitude of 
each mode is less than about 15 /im/s. It has been found 
that effects of " critical velocity" and hysteretic behavior 
are small at these low velocity amplitudes. The shifts 
in frequency and amplitudes of the modes can strongly 
affect each other if the drive levels are increased beyond 
some threshold velocityi-® 

Prior to loading the cell with a sample solid, the "back- 
ground" characteristics of the TO with no helium in 
the sample chamber are measured during slow (~ 15 
mK/hr) cooling and verified for reproducibility by fol- 
lowing the temperature dependence of the resonant fre- 
quencies, /ib(T), and the oscillation amplitudes, Aii,{T), 
of both modes (z = 1, 2). The background quality factor, 
Qib{T), is computed from the tracked oscillation ampli- 
tude via Q^biT) = iA,b{T)/AM))Q^b{Tb), where the 
reference quality factor Qib(Tb) is determined from mea- 
sured exponential ring down time (reproducible within 
10 %) at Tb « 200 mK. 

Background characteristics measured at T = 30 mK 
are: fu = 493 Hz, /zh = 1164 Hz, Qib ~ T.TxlO^ and 
Q2b ~ 4.7x 10^. The observed background frequencies are 
in agreement within ^12 % with those computed from 
estimated moments of inertia and torsion constant of the 
rods. 

After a sample of solid ^He is grown as described above, 
the resonant frequencies decrease owing to the added mo- 
ment of inertia. "Loading frequency" A/? in the "zero" 
temperature limit is defined as A/f = fib — fis corre- 
sponding to the decrease in the mode frequency from the 
background to sample-filled cell measured near 30 mK. 
The measured values are A/° = 0.40 Hz and and A/2 
— 1.20 Hz for the 2:3 = 0.3 ppm sample shown in Fig. [2] 
These values are within 10 % of those estimated from the 
added inertia of solid ^He in the sample chamber. The 
variation in Aff from sample to sample is less than 3 %. 



After a solid ^He sample is grown in the cell, the same 
procedure as in the background characterization includ- 
ing the drive levels is followed to measure fis{T), Ais{T), 
Qis(Tb) and Qis(T). The sample is warmed up at most to 
300 mK in the a::3 ~ 0.3 ppm sample, to 1.1 K in the 3, 6 
and 12 ppm samples, and to 2.0 K in the 25 ppm sample. 
Each sample data set is taken over about 24 hour pe- 
riod while cooling down to ~20 mK from 250 mK in the 
0.3 ppm sample and from 1 K in all other samples. The 
measured temperature dependence of fis (T) and Ais (T) 
does not vary significantly from one temperature sweep 
to another if the maximum temperature is within these 
limits. 



III. RESULTS 
A. frequency shift 

Our results on frequency shifts are presented in Fig. [2] 
as "reduced frequency shifts" defined for each mode i as: 

. l(fam-A/?)-/..(r)] _ 

</ is 

where is the frequency of loaded TO at our mini- 
mum temperature (about 15 mK) depending on X3. In 
the zero temperature limit the reduced frequency shift 
vanishes by definition. It remains, in all samples, at the 
zero temperature limit below about 40 mK and its mag- 
nitude monotonically increases at higher temperatures 
for both modes. Except in the X3 = 0.3 ppm sample, 
temperature dependence of the reduced frequency shift 
of the first mode coincides with that of the second mode 
below about 100 mK. At temperatures greater than 100 
mK changes in reduced frequency of the second mode are 
greater than those of the first mode. In all samples, ex- 
cept possibly the 0.3 ppm sample, the reduced frequency 
shifts for both modes continue to decrease as tempera- 
ture is increased above 200 mK. The lack of data in the 
0.3 ppm sample prevents us from making a firm state- 
ment about the temperature dependence of the reduced 
frequency shift above 200 mK in this sample. 

The temperature range 50 mK ^ T 150 mK where 
relatively rapid changes in the reduced frequency shift oc- 
curs in Fig. [2] has been identified as a signature of the oc- 
currence of NCRI phenomenon. The reduced frequency 
shifts would become constant if the observed temperature 
dependence of fis (T) matched that of fib (T) at temper- 
atures greater than some " onset" temperature where the 
fraction of solid sample apparently decoupled from the 
container, or NCRI fraction (NCRIf), vanishes. Since 
identifying such onset temperatures in our samples is am- 
biguous, reduced frequency shifts rather than NCRIf are 
shown in Fig. [5] 

The reduced frequency shift in the 25 ppm sample 
shows qualitatively distinct behavior from other samples. 
In comparison to other samples, there is no temperature 
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FIG. 2. (color online) Temperature dependence of reduced 
frequency shift of the compound torsional oscillator with an- 
nular sample chamber: first mode ((red) open circles) and 
second mode ((blue) pluses). Panels (a) ~ (e) show results of 
solid *He samples containing ^He impurity concentration of 
0.3, 3, 6, 12 and 25 ppm, respectively. Indicated solid pres- 
sures are estimated from the loading liquid pressure at 4.2 
K. 
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FIG. 3. (color online) Temperature dependence of the change 
in dissipation from empty to loaded sample chamber: first 
mode ((red) open circles) and second mode ((blue) pluses). 
Measurements are taken simultaneously with those reduced 
frequency shifts in the samples (with the same panel desig- 
nations for '^He impurity concentration) as shown in Fig. m 
Arrows indicate temperatures where peaks in dissipation oc- 
cur. 



range where the reduced frequency shifts vary relatively 
more rapidly. At temperatures roughly above 150 mK 
the frequency shifts decrease linearly (on the logarithmic 
temperature scale) with the second mode having larger 
slope than the first. Measured dissipation in the 25 ppm 
sample is also distinct from other samples (see Fig. [S]). 

B. dissipation 

The change in dissipation of each mode produced by 
loading solid ^He samples is computed by taking the dif- 
ference: 

AQ7\t)^Q7/(T)-Q;,\T). (2) 

Evaluated temperature dependence of AQ^^(T) is dis- 
played in Fig.[3]for each sample shown in Fig. [2] Samples 
with 2:3 < 12 ppm have similar temperature dependence 
with AQY^{T) passing through local maxima at "peak 
temperatures" (Tip) indicated by arrows in Fig. [31 The 
25 ppm sample shows broader temperature dependence 
than the other samples. The change in dissipation of the 
first mode is greater than that of the second in all sam- 
ples. This is in contrast to the reduced frequency shift of 



the second mode being greater than that of the first in 
all samples at temperatures above 100 mK (see Fig. [5]). 
At temperatures below 40 mK, AQ^^ for some samples 
is negative. This peculiar feature may have resulted in 
part from an uncertainty (~10 %) in the measurements 
of reference quality factors and from small inaccuracy in 
the cell temperature measurement below 40 mK. Small 
temperature gradients between the sample and the ther- 
mometer could lead to errors in taking the difference be- 
tween the two measurements with the sample chamber 
being empty and loaded. 

The observed frequency dependent peak temperature 
that the compound TO technique uniquely yields is an 
important parameter in considering the dynamical effects 
occurring in the oscillating ^He samples. The peak tem- 
peratures are fairly well-defined except in the 6 and 25 
ppm samples whose AQ^^ are broader than the others. 
The inverse of the peak temperature of the two modes 
in each sample is shown in Fig. ID It can be seen that 
Tp2 > Tpi in all samples. This frequency dependence 
of the peak temperature was also found in a cylindri- 
cal sample chamber geometry in solid *He with nominal 
0.3 ppm ^He concentration^" and appears to be indepen- 
dent of sample geometry and size (down to 0.2 mm), '^He 
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FIG. 4. (color online) Inverse of dissipation peak tempera- 
ture vs. ■^He impurity concentration: first mode ((red) circles) 
and second mode ((blue) squares). "Half-maximum temper- 
atures" or Tso, (crosses) are taken from Kim, et al.— 



concentration and sample growth conditions. The peak 
temperature generally increases for both modes as ^3 is 
increased. An apparent plateau in Tip in the range 3 
ppm < 2:3 < 10 ppm is likely caused by some accidental 
variation in as yet unidentified source of dissipation in 
the sample characteristics. It is interesting to note that 
the "half-maximum temperature" (Tso), where Kim, et 
alj^ found NCRIf to decrease to half of the maximum at 
lowest temperatures, smoothly extends the dependence 
of Tip on X3 found here. 

The same results as shown in Fig. [3] are replot- 
ted in Fig. [5] by normalizing the inverse temperature 
as T-'^/T-^ and the dissipation as l\Ql'^ / /\Q,{T,p)~^ . 
Plotting in this manner reveals similarity and dissimilar- 
ity among the samples with varying ^He impurity concen- 
tration. The normalized dissipation for the first mode in 
the 12 and 25 ppm samples and the second mode in the 
25 ppm sample deviate considerably from others. The 
"sharp" increases in the normalized dissipation at high 
temperatures, Tip/T ^ 0.2 (first mode) and 0.4 (second 
mode), become accentuated in Fig. [5l Deviations in 
the samples with higher X3 from the other bell-shaped 
dependence on Tip/T indicate an emergence of nearly 
temperature independent "extra" dissipation. 

Annealing of solid ^He samples has been found^i^i^ 
to affect the frequency shift and dissipation in TO ex- 
periments. Possible effects of sample "annealing" were 
studied in our X3 — 0.3 ppm sample by raising the tem- 
perature of the cell up to 1.8 K (below the melting tem- 
perature of the sample) and maintaining the temperature 
for 10 hours. The temperature is then reduced to 0.3 K 
over six hours. Subsequently measured temperature de- 
pendence of the frequency shifts of the two modes are 
similar to those shown in Fig. [2] except /£, decreases by 
7.2 and 21.5 mHz for j = 1 and 2, respectively. This ob- 
servation is contrary to that of the recent experiment^^ 
which finds that the measured frequency shift at lowest 
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FIG. 5. (color online) Normalized (see text) change in dissipa- 
tion: first mode (panel (a)) and second mode (panel (b)) of all 
samples shown in Fig. [3l Symbols: (purple) circles (3:3 [ppm] 
= 0.3), (red) lozenges (3), (green) triangles (6), (blue) squares 
(12), and (orange) stars (25). Results shown for the 25 ppm 
sample include temperatures up to 2.5 K not displayed in Fig. 
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temperatures is not affected by annealing. The dissipa- 
tion peak temperatures decrease by about 5 mK after an- 
nealing. Our annealing process decreases the extracted 
(see Sec. lIVp characteristic time slightly but produces no 
significant change in the activation energy. 



IV. ANALYSIS 

The original discovery- of frequency shifts at low tem- 
peratures in TOs loaded with solid ^He and all subse- 
quent confirmations, to our knowledge, have been ac- 
companied by dissipations having puzzling resonance- 
like temperature dependence as exemplified by panel 
(a) in Fig. [31 In this section, we analyze the ob- 
served dependence of dissipations and frequency shifts 
on ^He impurity concentration in solid '^He samples. 
It is generally agreed that dissipation peaks occur at 
temperatures, where an internal dynamical rate (t~^) 
matches the imposed TO frequencies (wi), or where 
ujiT = 1. The internal dynamics of vortex motion!^, 
glassy response^, superglass^i^, presence of tunneling 
two level level systems^, and viscoelastic behvior— have 
been suggested as the physical origin of dissipation. The 
vibration of dislocation line segments pinned at network 
nodes and by '^He impurity has also been suggested^ as 
the origin of the observed TO behavior. 
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FIG. 6. (color online) Torsional oscillator frequencies plot- 
ted at inverse temperatures where dissipation peaks occur in 
solid "^He samples with given "^He impurity concentrations. All 
straight hues have the same slope corresponding to an activa- 
tion energy Eo of 430 mK and have ordinate intercepts which 
are adjusted for best fit for each sample. The intercepts repre- 
sent the characteristic dynamic rate rj~^ (see text). Symbols 
for samples with different values of X3 are same as those in 
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FIG. 7. Characteristic time tq vs He impurity concentra- 
tion. The values of tq are determined by the intercepts of 
the straight lines in Fig. [6l The dashed line represents 
To — 3.5 X 10"''' 2:3''^ s (see Discussion section of text). Crosses 
show To calculated from the T50 data of Kim, et alp^ using the 
same activation energy as in Fig. [6] and their TO frequencies 
(see text). 



Assuming that the internal dynamics is thermally 
driven, we consider a simple, activated dynamical time t 
given by the Arrhenius form: 



T = Tq exp 



En 



ksiT-ToY 



(3) 



where fcs is the Boltzmann constant, tq a characteris- 
tic time and Eq an activation energy. A possibly non- 
vanishing "transition temperature" To is introduced in 
Eq. ^ for generality, but we assume Tq = in our 
analysis for simplicity. To examine our dissipation re- 
sults in terms of Eq. ([3|), the mode frequencies are 
plotted in Fig. |4] on a logarithmic scale at correspond- 
ing in each sample. The slopes and the intercepts 
determined strictly by straight lines {not shown) con- 
necting the two points for each sample in Fig. [S] are 
(i;o/fcs[mK],To-i[s"i]) = (380, 5x10"^), (290, TxlG-^), 
(720, 2x10-8), (370, SxlO^^) and (430, 3x10-^) for xz 
= 0.3, 3, 6, 12, 25 ppm samples, respectively. Broadness 
in the dissipation peaks introduces considerable uncer- 
tainties (indicated by error bars in Fig. in the values 
of Tip and leads to scattering in these values of Eo/ks 
and Tq^. The tendency in Fig. |6l however, suggests that 
Eq is a constant independent of X3. To make progress, let 
us assume that Eq/Ub = 430 mK, the average of above 
slopes, is a good estimate for the activation energy. It is 
interesting to note that this average activation energy is 
close to the binding energy of '^He dislocation lines found 
in the analysis of their experiments by Kim, et al.— and 
Day and Beamish.— 

The intercepts of the best fits (shown by a straight line 
for each X3 sample in Fig. [5]) with the set slope specify 



the values of tq as shown in Fig. [Tj There is a trend for 
Tq to increase as 0:3 where 7 ~ 2/3 if X3 < 20 ppm. Both 
the reduced frequency shifts and dissipation of 25 ppm 
sample show qualitatively different behaviors than the 
samples with lower '^He impurity concentration. It ap- 
pears that Tq also changes its dependence on beyond 
20 ~ 25 ppm. It was already noted that the values of T^q 
measured by single frequency TO techniques by Kim, et 
al.Si are fairly close to Tip (see Fig. 3]). Constraining 
straight lines with the same slope as shown in Fig. [Hby 
their TO frequencies at T^^ imply values of tq as shown 
in Fig. [7] for their samples. Considering differences in 
sample chamber geometry, sample growth process, mea- 
surement methods, etc., tq extracted from T^q data of 
Kim, et al.?i overlaps and fits surprisingly well with ex- 
trapolation of our results to lower values of x^ . Clearly 
it is of interest to extend frequency dependence studies 
like ours to smaller "^He impurity concentration than we 
have carried out. 

The mechanical response of driven TOs containing 
solid ^He samples has been treated^! by including a gen- 
eral form of rotational susceptibility to account for the 
sample motion: 



2G 



1 - [jUJiT) 



(4) 



where G is a constant (possibly dependent on 
frequency^), uji is angular frequency, r is a relaxation 
time, and j3 is an exponent dependent on a particular 
model. Including the general susceptibility results in 
extra dissipations and concomitant frequency shifts in 
the TO response. In our simplified analysis, the "Debye 
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model" response is considered by assuming /3 = 1. In 
this case, the change in dissipation due to sample motion 
is given by: 



(l+c.fr2)' 



and the accompanying reduced frequency shift by: 



fis fib ^fi 



G 



(1 



(5) 



(6) 



Although the linear response model assumed here may 
not be entirely applicable to our experiment, it is of inter- 
est to examine if the model can achieve a similar success 
as in the description^ of the shear modulus measure- 
ments. Our results are compared with those expected 
from Eq. ([5]) by assuming the Arrhenius form of relax- 
ation time characterized by our average activation energy 
and the extracted tq as shown in Fig. [T] Dissipations 
extracted from the experiment and from the model are 
compared in Fig. |5]for the 0.3 ppm sample. The observed 
peak dissipation of the first mode at Tip is significantly 
larger than that of the second mode at T2p in disagree- 
ment with the frequency independent peak dissipation 
expected from Eq. ([5]). To proceed with the (modified) 
Debye model, the value of G is adjusted for each mode 
separately to match the peak dissipation at Tip in Fig. |51 
Dissipations expected from the model produce narrower 
widths in temperature dependence around the peaks than 
those observed. Comparisons in other samples show sim- 
ilar deviations between the observations and the model. 
The model of course does not account for the upturns 
in dissipation at high temperatures nor for the low tem- 
perature residual dissipations observed in 12 and 25 ppm 
samples. 

As seen in Fig. [SI the observed temperature depen- 
dence of dissipation is broader than than expected from 
the Debye model assuming Arrhenius relaxation time 
with one activation energy for the system. Broader tem- 
perature dependence may be introduced into the model 
by allowing a distribution in the activation energy. A 
canonical Gaussian distribution N{E) given by 



N{E) = 



1 



1 , E-Eq ^2 



is applied to Eq. ^ to evaluate AQ~^: 

AQr\T)= [ AQr\T,E)N{E)dE. 



(7) 



(8) 



Here, Eo/ks is set as the average activation energy and 
the width w of the distribution is adjusted for each sam- 
ple. The peak dissipation value G is readjusted separately 
for each mode to match the respective peak dissipation. 
For the 3 ppm sample, adjusting to w/ks = 120 mK can 
represent the data fairly well as shown in Fig. [Sj In 
the cases of 12 and 25 ppm samples, there appear resid- 
ual amounts of dissipation at low temperature. In these 
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1000 



FIG. 8. (Color online) Comparison of dissipation in the 3 ppm 
sample with the dissipation expected from Debye susceptibil- 
ity. Same data as shown in Fig. [3] are repeated for mode 
1 ((red) circles) and 2 ((blue) pluses) in panel (a) and (b), 
respectively. In each panel, curves are dissipations expected 
(see text) from Eq. ((5|: assuming single activation energy and 
To = 1.2xl0~® s and from Eq. ((Tj) (dashed curves), and al- 
lowing Gaussian distribution of activation energy with w/ks 
= 120 mK and from Eq. ((8| (solid curves). 



250 



200 



150 




100 



(ppm) 

FIG. 9. Variation of width in distribution of activation energy. 
Widths are determined by fitting temperature dependence of 
dissipation in each sample using Eq. (O and ((§]) (see text). 



two samples, a constant added to AQ~^ is taken as an 
additional fitting parameter. The width of distribution 
obtained from fitting the data in this manner is shown in 
Fig. ini Values of w within the error bars shown in Fig. 
ini give similar goodness of fit for both modes. Despite 
large uncertainties, a clear increasing trend in w is dis- 
cernable as is increased. We suggest that temperature 
widths around the dissipation peaks in TO experiments 
can provide sample characterization in the distribution 
of activation energy. 

Using the same respective parameter values in eval- 
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uating dissipation for each mode shown in Fig. |8l the 
frequency shifts expected from Eq. with single and 
distributed activation energy, are shown in Fig. [inifor the 
3 ppm sample. The expected frequency shift of the first 
mode by including distribution in the activation energy 
approaches the observed frequency shift but a consider- 
able difference remains. The expected frequency shifts 
of the second mode are much smaller than the observed. 
The large discrepancy in the second mode occurs despite 
independently adjusting the value of G for this mode as 
described above. The observed decrease in dissipation 
and increase in reduced frequency shift as the TO fre- 
quency is increased is a major inconsistency with the 
modified Debye model above. The difference in reduced 
frequency shifts between the measured and the Debye 
model would be a consequence of a superfluidity in solid 
^He. A similar conclusion is made by Yoo and Dorsey^ 
in the analysis of their viscoelastic model for solid ''He. 

Despite the clear distinction demonstrated in the re- 
duced frequency shifts between the measured and the 
modified Debye model, similarity in temperature depen- 
dence between the two in Fig. [10] is evident. Multi- 
plying the Debye expectations shown by solid curves for 
the first and second mode by constants, 1.7 and 6.0, re- 
spectively, gives temperature dependence shown by dash- 
dotted curves in Fig. [TUl Except in high temperature 
tail regions, the scaled temperature dependence matches 
the measured reduced frequency shifts quite well. Similar 
matching is seen in other samples with multiplicative con- 
stants applied to (first mode, second mode) of (1.4, 3.5), 
(1.3, 3.9), and (1.25, 3.6) in 0.3, 6 and 12 ppm samples, 
respectively. The 25 ppm sample cannot be matched in 
a similar manner owing to the large offset in dissipation. 
How this similarity in temperature dependence is related 
to the interpretation of supersolidity is not yet clear to 



Nussinov, et alM- and Graf, et alM have carried out 
much more sophisticated analyses of observed TO re- 
sponses with generalized rotational susceptibilities ex- 
pressing glassy response. They found that various sin- 
gle frequency TO responses in both dissipation and fre- 
quency shift could be fitted by appropriately adjusting 
parameters in the glassy response model. Graf, et ali^ 
recently reported their study of the frequency depen- 
dence in the dissipation and frequency shifts observedi^ 
in our compound TO with a cylindrical sample chamber 
by allowing frequency dependence in the parameter G in 
Eq. (j?]). In our simple analysis of the observed ^ He impu- 
rity concentration dependent effects in the present annu- 
lar sample chamber, we restricted ourselves to the above 
modified Debye model in which the value of G is adjusted 
for each mode separately to account for its apparent fre- 
quency dependence. Gomparison of the cylindrical and 
annular sample chamber results indicates that the magni- 
tude of frequency dependence of reduced frequency shifts 
increases as the sample chamber size is decreased. 




100 
T(mK) 

FIG. 10. (color online) Comparison of 2:3 = 3 ppm sample 
frequency shift data with those expected from Debye suscep- 
tibility using parameters determined from fitting the dissipa- 
tion data as shown in Fig. [S] Dashed and solid curves are 
those expected from the dissipations shown in Fig. [8]with the 
single activation energy and with the distribution of activa- 
tion energy, respectively. Dash-dotted curves are obtained by 
multiplying the solid curves of the respective modes by con- 
stants (see text). Same data (but depopulated for clarity) as 
shown in Fig. [2] are repeated for mode 1 ((red) circles) and 2 
((blue) pluses) in panel (a) and (b), respectively. 



V. DISCUSSION 

It is noted that the extracted characteristic time shown 
in Fig. [7] varies almost in proportion to x^^^ except for 
the 25 ppm sample. We speculate that this dependence 
stems from the diffusion process of '^He condensed onto 
dislocation lines in solid ^He samples. The diffusion time 
Td of •^He along the dislocation lines may be approxi- 
mated as Td « /D, where s is some characteristic dis- 
tance over which "^He moves during a time interval 
and D K Iv is a, diffusion constant with mean free path 
I and particle velocity v. Let us suppose that the mean 
free path is approximately given by the effective disloca- 
tion loop length L. We follow Iwasa^ in writing L as a 
parallel combination of the "^He impurity length Li and 
the network pinning length Ljv: 



L 



N 



gx^ ^ exp 



2Wo 
' 3T 



L 



N ' 



(9) 
(10) 



where g — 3.4 x 10~^(see Iwasap^^) is a constant and Wq 
is the binding energy of '^He impurity to dislocation line. 
At low temperatures where Li dominates over Ljv, the 
characteristic time takes on the Arrhenius form: 



Td = — x§ exp(£'o/T) 
gv 



(11) 



Identifying Td as r and the coefficient on the right hand 
side of Eq. PT|) as tq, we expect it to depend on the 
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impurity concentration as cx . Letting s^/u = 1.2 x 
10-9 cm s, Wo/fcs (= 3Eo/2kB) = 0.65 K, = 2 ^m 
(similar "fits" can be achieved for 1 ^ 5 /im) gives the 
dependence on 2:3 as shown by the dashed hne in Fig. 
[71 This dependence appears to describe the observed 
dependence of tq on X3 except for the 25 ppm sample. 

There are questions that can be raised on the above 
diffusion process to describe the observed characteris- 
tic time To as X3 is varied. If the mean free path of a 
•^He atom along a dislocation line is determined by other 
•^He atoms condensed onto the dislocation line, it implies 
that ''He atoms along the line mysteriously manage not 
to contribute to scattering. The process also implies that 
the diffusion coefhcient of "^He would decreases in the low 
temperature limit as exp (—Eq/T). This temperature de- 
pendence is in contrast to the temperature independent 
quantum diffusion coefficient of '^He in bulk solid ''He 
found by NMR experiments^'^ above 0.55 K. The Kyoto 
group has observed very long spin-lattice relaxation in 
their NMR experiment on solid ''He sample with 2:3 w 30 
ppm. The long relaxation time might be related to the 
decreasing D implied by our analysis. 

Iwasa— has analyzed TO experiments and shear mod- 
ulus shifts in terms of the Granato-Liicke model^ on 
the interaction between an externally oscillated TO con- 
tainer and the induced vibrational motion of the dislo- 
cation lines present in the loaded ^He solid sample. The 
analysis predicts shifting temperature dependence of the 
TO frequency and the dissipation peak temperature to 
higher temperatures as X3 is increased in general agree- 
ment with observations. However, since the expected 
natural vibration frequencies of the dislocation lines are 
much higher than the TO frequencies so far attempted, 
little frequency dependence is expected in both frequency 
shift and dissipation. This aspect of the model^^ is yet 
to be reconciled with our TO experiments. 

Gaudio, et al4^ proposed a model*^ on the effects 
of ■^He on TO experiments where uniformly distributed 
■^He impurities set the maximum grain size in solid "^He 
samples. This model is not likely to be applicable in 
the range of values of 0:3 in our experiments where '^He 
atoms are expected to condense onto dislocation lines. 
Manousaki^ considered '^He impurity atoms binding to 
defects and promoting "^He atoms to interstitial sites. It 
is not clear to us how this effect relates to the temper- 
ature dependent changes in dissipation observed here as 
2:3 is varied. 

Day and Beamish^'^ discovered that temperature de- 
pendent changes in the shear modulus of solid "^He were 
almost identical to those in frequency shifts observed 
in TO experiments. The discovery gives strong impe- 
tus for concluding that the observed changes in shear 
modulus and frequency shift have similar physical origin. 
Syshchenko, et al.-^ recently reported on their measure- 
ments of the changes in shear modulus and associated 
dissipation at frequencies between 0.5 Hz and 8 kHz. 



They analyze (see also Su, et al.-^) their data with a De- 
bye relaxation process and thermally activated dynamics 
just as we do. When a distribution of activation energy 
is included, they find that the relationship between the 
shear modulus and the dissipation can be accounted for 
at different frequencies by the assumed Debye and ther- 
mal activation processes. This is in sharp contrast to the 
analysis of our TO data: including distributions in ac- 
tivation energy cannot account for both frequency shifts 
and dissipations. The discrepancy increases at higher 
frequency (see Fig. [T0|) . 

VI. CONCLUSION 

Effects of adding ■^He impurity to solid '*He sam- 
ples contained in an annular chamber were studied si- 
multaneously at two resonant mode frequencies (differ- 
ing by a factor of 2.4) by means of a compound tor- 
sional oscillator. Both frequency shifts and extra dis- 
sipations produced by the loaded samples were mea- 
sured. Maxima in the measured dissipation occurred 
at impurity-concentration (0:3) dependent "peak temper- 
atures" around which the frequency shifts varied more 
rapidly. When normalized to both the temperature and 
the dissipation level at the peak, the temperature depen- 
dence of dissipation became nearly universal in all sam- 
ples studied except in the 25 ppm sample. A thermal 
activation energy (430 mK) and characteristic relaxation 
times were extracted from Arrhenius plots of frequency 
versus the inverse of dissipation peak temperature. The 
characteristic time increased with impurity concentration 
approximately as x^^^ and suggested diffusion of '^He 
atoms along dislocation lines as the dynamical process 
producing the observed dissipation. Observed tempera- 
ture dependence of dissipation of both modes could be 
fairly well described by a simple Debye model by allow- 
ing for Gaussian distribution of activation energy if the 
magnitude of dissipation at the peak temperature was 
allowed to be frequency dependent. The measured mag- 
nitudes of frequency shifts were significantly greater than 
those expected from the model especially in the higher 
frequency mode. There remained "excess" amounts of 
frequency shifts which could not be accounted for by the 
simple Debye model. The excess frequency shifts may be 
attributed to superfluidity in solid ^He at low tempera- 
tures. We believe that these frequency dependent effects 
hold a key to understanding of the dynamics of quantum 
solid ^He. 



VII. ACKNOWLEDGMENT 

We thank Moses Chan for stimulating discussions. 
This research was supported by NSF DMR-0704120 and 
DMR-1005325. 



10 



* Present address: Institute for Quantum Computing, Uni- 
versity of Waterloo, ON, Canada. 

1 E. Kim and M. H. W. Chan, Nature(London), 427, 225 
(2004); Science, 305, 1941 (2004). 

2 N. Prokof'cv, Adv. Phys., 56, 381 (2007). 

^ S. Balibar and F. Caupin, J. Phys. Cond. Mat., 20, 173201 
(2008). 

* D. E. Galli and L. Reatto, Journal of the Physical Society 
of Japan, 77, 1010 (2008). 

^ A. J. Leggett, Phys. Rev. Lett., 25, 1543 (1970). 
^ M. W. Meisel, Physica B, 178, 121 (1992). 

A. S. C. Rittncr and J. D. Reppy, Phys. Rev. Lett., 97, 

165301 (2006). 

* M. Kondo, S. Takada, Y. Shibayama, and K. Shirahama, 
J. Low Temp. Phys., 148, 695 (2007). 

® A. Penzev, Y. Yasuta, and M. Kubota, J. Low Temp. 

Phys., 148, 677 (2007). 
^° Y. Aoki, J. C. Graves, and H. Kojima, Phys. Rev. Lett., 

99, 015301 (2007). 
" B. Hunt, E. Pratt, V. Gadagkar, M. Yamashita, A. V. 

Balatsky, and J. C. Davis, Science, 324, 632 (2009). 

D. S. Greywall, Phys. Rev. B, 16, 1291 (1977). 

J. Day and J. Beamish, Nature(London), 450, 853 (2007). 

A. S. C. Rittncr, W. Choi, E. J. Mueller, and J. D. Reppy, 

Phys. Rev. B, 80, 224516 (2009). 

M. W. Ray and R. B. Hallock, Phys. Rev. Lett., 100, 
235301 (2008). 

^® M. W. Ray and R. B. Hallock, Phys. Rev. Lett., 105, 
145301 (2010). 

A. F. Andreev and I. M. Lifshitz, Sov. Phys. JETP, 29, 
1107 (1969). 

^® M. R. Sears and W. M. Saslow, Phys. Rev. B, 82, 094519 
(2010). 

Y. Aoki, H. Kojima, and X. Lin, Low Temp. Phys., 34, 
329 (2008). 

^° S. Kwon, N. Mulders, and E. Kim, J. Low Temp. Phys., 
158, 590 (2010). 

21 E. Kim, J. S. Xia, J. T. West, X. Lin, A. C. Clark, and 
M. H. W. Chan, Phys. Rev. Lett., 100, 065301 (2008). 

22 R. Toda, P. Gumann, K. Kosaka, M. Kanemoto, W. Onoe, 
and Y. Sasaki, Phys. Rev. B, 81, 214515 (2010). 

23 R. Toda, W. Onoe, M. Kanemoto, T. Kakuda, Y. Tanaka, 
and Y. Sasaki, J. Low Temp. Phys., 162, 476 (2011). 



2* S. S. Kim, C. Huan, L. Yin, J. Xia, D. Candcla, and N. S. 

Sullivan, J. Low Temp. Phys., 158, 584 (2010). 
25 S. S. Kim, C. Huan, L. Yin, J. Xia, D. Candela, and N. S. 

Sullivan, ArXiv e-prints/10123709. 
2^ Y. Aoki, J. Graves, and H. Kojima, J. Low Temp. Phys., 

150, 252 (2008). 
2^^ Emerson and Cuming, Billcrica, MA 01821. 
2* M. C. Keiderling, P. Gumann, D. Ruffner, and H. Kojima, 

To be published. 
2'' J. D. Reppy, Phys. Rev. Lett., 104, 255301 (2010). 
3° P. W. Anderson, Nat. Phys., 3, 160 (2007); see also S. L 

Shcvchcnko, Sov. J. Low Temp. Phys., 13, 61 (1987). 

31 Z. Nussinov, A. V. Balatsky, M. J. Graf, and S. A. Trug- 
man, Phys. Rev. B, 76, 014530 (2007). 

32 G. BiroU, C. Chamon, and F. Zamponi, Phys. Rev. B, 78, 
224306 (2008). 

33 A. F. Andreev, JETP, 109, 103 (2009). 

34 C. D. Yoo and A. T. Dorsey, Phys. Rev. B, 79, 100504 

(2009) . 

35 L Iwasa, Phys. Rev. B, 81, 104527 (2010). 

3® M. Graf, Z. Nussinov, and A. Balatsky, J. Low Temp. 

Phys., 158, 550 (2010). 
3'' O. Syshchenko, J. Day, and J. Beamish, Phys. Rev. Lett., 

104, 195301 (2010). 

3® M. J. Graf, A. V. Balatsky, Z. Nussinov, I. Grigorenko, and 
S. A. Trugman, J.Phys. Conf. Scries, 150, 032025 (2009). 

3^^ M. Graf, J.-J. Su, H. Dahal, I. Grigorenko, and Z. Nussi- 
nov, J. Low Temp. Phys., 162, 500 (2011). 

'*° 1. Iwasa and H. Suzuki, J. Phys. Soc. Jpn, 49, 1722 (1980). 

''1 A. R. Allen, M. G. Richards, and J. Schratter, J. Low 
Temp. Phys., 47, 289 (1982). 

''2 A. Granato and K. Liickc, J. App. Phys., 27, 583 (1956). 

43 S. Gaudio, E. CajipoUuti, G. RastcUi, and L. Pietronero, 
Phys. Rev. Lett., 101, 075301 (2008). 

4"* B. Yucesoy, J. Machta, N. Prokofev, and B. Svistunov, 
Phys. Rev. Lett., 104, 049601 (2010); S. Gaudio, E. Cap- 
pelluti, G. Rastelli, and L. Pietronero, ihid., 104, 049602 

(2010) . 

45 E. Manousakis, Europhysics Letters, 78, 36002 (2007). 
4« J.-J. Su, M. J. Graf, and A. V. Balatsky, Phys. Rev. Lett., 

105, 045302 (2010). 



